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Abstract 

The calculation of the 0(Nfa 2 ) corrections to muon decay 
is described. These are the 2-loop diagrams containing a mass¬ 
less fermion loop and they form an important gauge-invariant 
subclass. It is shown that all such diagrams can be expressed 
in terms of a universal master integral. We focus on the cal¬ 
culation of box diagrams and in particular on the removal of 
their infrared divergences. 


1 Introduction 

Muon decay, /i~ —> u^e~v e , has always been a proving ground for 
both pure QED and electroweak radiative corrections [jTj-0. Until re¬ 
cently the muon decay coupling constant, G^, and the electromagnetic 
coupling constant, a , were by far the best measured electroweak ob¬ 
servables and played a pivotal role as input to the Standard Model. 
The accuracy of theoretical predictions was limited by the errors on 
Mz which was taken as the third input required to make the model 
predictive. Now the situation has changed somewhat. Both Mz and 
G M are determined to an accuracy of 2 x 10~ 5 and it may be possi¬ 
ble to reduce the error on Mz still further if the LEP beam energy 
calibration can be better understood. 

By comparison electroweak calculations have attained an preci¬ 
sion of only a few permill. Obviously it would be desirable to have 
the results of full 2-loop electroweak calculations at our disposal and 
considerable progress has been made in this direction JiJ. Such cal¬ 
culations are likely to be intractable when expressed in analytic form 
and in practice one would quote only numerical results. There do ex¬ 
ist, however, a few 2-loop calculations for which exact analytic results 
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can be obtained in a compact closed form IBM Here we discuss 
one such set of corrections. These are dubbed 0(Nfa 2 ) corrections to 
muon decay, where Nf is the number of light fermions. Because Nf is 
quite large, the 0(Nfa 2 ) corrections are expected to be a dominant 
subset of 2 -loop graphs and since Nf provides a unique tag, the com¬ 
plete set of (D(Nfa 2 ) corrections contributing to a particular physical 
process will form a gauge-invariant set. Diagrams of this type have 
been considered in connection with the muon anomalous magnetic 
moment 0 . 

A priori the 0(NfOi 2 ) corrections can be expected to contribute 
at the level of 1.5 x 10~ 4 . Without their calculation and inclusion, 
theoretical predictions remain uncertain at least at this level. 

If the fermions are assumed to be massless compared to Mw, an 
excellent approximation for all but the third generation, then the only 
‘dimensionfuh quantities appearing in the calculation are Mw , Mz 
and Mh and there are very few diagrams containing a Higgs. The 
contributions from diagrams not involving Higgs bosons can then only 
be a polynomial in sin 2 6w with coefficients involving lncos 2 0 vv; hi 7 r 
and Euler’s constant, 7 . 

Assuming massless fermions reduces the topologies of diagrams 
that must be considered since the fermions can then only couple to 
vector bosons but not to Goldstones. However, the calculation is still 
a fully-fledged 2 -loop electroweak calculation requiring the complete 
renormalization of the Standard Model at 0(Nfa 2 ). Although muon 
decay represents a zero momentum transfer process, much can be 
learned from it about the cancellation of divergences in high-energy 
processes. Indeed unlike the case of the calculation of 0(a 2 m 4 ) cor¬ 
rections to the p-parameter ||j(|.| 8 |], there is a proliferation of diagrams 
involving counterterms and these constitute significant fraction of the 
effort involved. 

In addition to the calculation of diagrams at zero momentum trans¬ 
fer one also requires the W and Z° mass counterterms that must be 
obtained by evaluating diagrams at high scales. In principle this can 
be done in any renormalization scheme using methods expounded in 
ref.s [| 10 |,ITIJ but the work is considerably reduced by adopting the MS 


scheme in which only the simpler divergent parts of the diagrams are 
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required. 

In this talk we will summarize some of the salient features of the 
calculation. In section 2 the master integral will be given and in sec¬ 
tion 3 we will concentrate on contributions coming from box diagrams. 
Diagrams of this type do not occur in other 2-loop electroweak calcu¬ 
lations performed to date. It will be shown how the IR divergences 
can be separated from these diagrams in a well-defined manner that 
permits them to be meshed with bremsstrahlung diagrams. 

The work reported here was undertaken with R. Akhoury and 
P. Malde. Full details of the general methods used can be found in 
ref. [12] and results for muon decay in ref. |13| . 


2 The Master Integral 


In what follows an anticommuting y 5 is assumed. A Euclidean metric 
with the square of time-like momenta being negative will be used and 
all calculations will be done in R^=i gauge. The sine and cosine of the 
weak mixing angle, 9w, will be denoted sg and eg respectively, and 
7 r are the usual left- and right-handed helicity projection operators. 

Many, but by no means all, 0(Nfa 2 ) diagrams are obtained simply 
by inserting a fermion loop into the boson propagator of a one-loop 
diagram. Because the original one-loop diagram is often logarithmi¬ 
cally divergent, the fermion loop insertion will need to be calculated 
to 0(n — 4) in dimensional regularization, where n is the dimension 
of space-time. It turns out, however, that for all 0(NfOe 2 ) diagrams 
occurring in low-energy processes, it is possible to obtain expressions 
that are exact in n. 

For the massless fermion loop insertion, it may be shown that 



(71 
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where (3l, /3r and /3' L , /3' R of the couplings of the attached vector 
bosons. 
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For processes at zero momentum transfer it is possible to immedi¬ 
ately reduce all tensor integrals that occur to scalar integrals by using 
projection operator techniques ||T 2 |. These scalar integrals can then 
be written as expressions involving a general master integral, 


d n p 


d n q 


in 2 [p 2 ] J '[p 2 + M 2 ] k J in 2 [g 2 ] ; [(g + p) 2 } 
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3 Box Diagrams 

A useful set of identities for calculating one-loop box diagrams appears 
in ref. [14]. They are, however, valid only for n = 4 because of their 


intended use at one-loop. For general n it may be shown that these 
relations become 


hwvyL,a] 

1 [lplvlalL,B\ 2 

= 

[7a7i,fl] 1 [7«7I/,r] 2 




+ 

(n — / Y)[1 p.1 L ,r\i[1v1 L ,r\i 

(3) 


l\lpi ulalR^ 

= 

4 [1v1l,r] 1 [1p1r,l\ 2 




+ 

(n ~ 4)[1p1l,r]i[1v1r,l] 2 

(4) 

[7p7/u7o-7l,r] 

1 [1o1v1 p 1l,r\ 2 

= 

^{ivl L,r\ lb P-l L,r\ 2 




+ 

(n -4) [1p1l,r] r [7^7 l,r] 2 

(5) 

[7p7m7<t7l,.r] 

1 [ivlvlplRl]? 

= 

4^1/ [1<x1L,r} 1 [lalRll 2 




+ 

(n ~ 4:)bp,l L MlulR, L ) 2 

(6) 


where the square brackets [ ]i and [ ] 2 indicate that the enclosed 7 - 
matrices are associated with the separate external fermion currents J\ 
and J 2 respectively. These identities, along with projection operator 
techniques, may be used to reduce box diagrams to the form I ■ A4o 
where I is a scalar integral and A4q is the Born level matrix element. 

All 0(Nj a 2 ) box diagrams correspond to the insertion of a fermion 
loop into boson propagators in one-loop box diagrams. The box dia- 
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Figure 1: A class of box diagrams containing UV and IR divergences. 


grains, so obtained, are either logarithmically divergent or have dou¬ 
ble poles at n = 4 corresponding to mixed UV and IR divergences or 
fermion mass singularities. 

Consider the diagrams of Fig.l. The crosses represent fermionic 
parts of 0(a ) counterterms. The fermion loop contribution in Fig.la 
vanishes at q 2 = 0 where q is the 4-momentum in the photon prop¬ 
agator and the diagram is thus IR finite. If the counterterms on the 
W line in Fig.lb-d together vanish at q 2 = 0 then they too will be 
IR finite however this depends on the renormalization scheme that is 
chosen and is not generally the case. It happens, for example, in the 
MS bnt not in the on-shell renormalization scheme. 

In Fig.lb-d the insertion of counterterms corresponds to replacing 
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Figure 2: A class of box diagrams containing UV and IR divergences 
and fermion mass singularities. 


the W propagator in the one-loop box diagram 


q 2 + 

(Jj\ t _ Ml f / Sg\ _ SMI 1 

V 9 ) (9 2 + miy ( 9 = + Mlf H g) Ml f 

( 8 ) 

where 5g and are the SU(2 ) coupling constant counterterm and 
W mass counterterm respectively. 

The first term of (§) yields an IR finite contribution and its UV 
divergences cancel against those of Fig. la. The second term in (^]) 
is UV finite but generates an IR divergence and must be combined 
with soft bremsstrahlung correction to produce a finite result. As 
mentioned above the last term vanishes in the MS scheme. 

A similar thing happens for the diagrams shown in Fig.2 except 
that now Fig.2a is both IR divergent and has a fermion mass singular¬ 
ity. If we denote the fermion loop insertion in the photon propagator 
as (g 2 <5^ - g /i g^)hh f7 (g 2 ) then we may write 

n; 7 (g 2 ) = [n; 7 (o)] + [n; 7 (o) - n; 7 (o)] + [n; 7 ( g 2 ) - n; 7 (o)] 

(9) 


5g\ 1 

g)q 2 + Ml {q 2 + M^Y 
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where 


fp /'n't _ 2^9 . 2^9 / in \ 

n 77 (o) — s e — + c e~~r (io) 

and g' is the U( 1) coupling constant. 

The hrst term on the rhs of (|9|) yields an IR divergent correction 
that cancels against the IR divergences of the counterterm diagrams 
Fig.2b and 2c leaving an overall UV divergence. This remaining UV 
divergence cancels against those of other box diagrams involving Z° 
bosons. The second term in square brackets yields an IR divergent 
term that needs to be combined with soft bremsstrahlung diagrams 
to yield a finite result. Finally the third term gives a finite contribu¬ 
tion that is singular for vanishing fermion masses. It will be subject 
to a hadronic contribution and should therefore be evaluated using 
dispersion relations. We do not consider the second and third term 
further. 

Once IR divergences are removed in this way, the 0(Nfa 2 ) correc¬ 
tions to the Born level matrix element, Ado, coming from box diagrams 
may be written as AdoAwhere 

Ar<2) = (lfj 3 ) + 12s "“ 131s » + 135s " ” 45 ' 

-hn4(2 si - 1 ) - ss[ng(-M|) - n; 7 (o)] 

2 1 (11) 
+h#[n g(-M|) - ng(o)](2sJ -104 + 5 ) 1 

s e J 

for each massless fermion generation. 

This result is given in the on-shell renormalization scheme and 
the photon vacuum polarization, IT (—M§) —111^ (0), arises from the 
definition of one-loop counterterms in this scheme. The correspond¬ 
ing expression in the MS scheme is somewhat longer because of the 
presence of terms such as In M^jg 2 and In 7r. 
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